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Abstract 

We investigate the ground-state properties of a collection of N non-interacting electrons in a 
macroscopic volume also containing a crystalline array of N spheres of radius rc each taken as 
largely impenetrable to electrons and with proximity of neighboring excluding regions playing a 
key physical role. The sole parameter of this quantum system is the ratio r^/rg., where is the 
Wigner- Seitz radius. Two lattices (FCC and BCC) are selected to illustrate the behavior of the 
system as a function of rdrg. As this ratio increases valence electrons localize in the interstitial 
regions and the relative band- width ep/e^ is found to decrease monotonically for both. The system 
is motivated by the behavior of the alkali metals at significant compression. It accounts for band 
narrowing, leads to electronic densities with interstitially centered maxima, and can be taken as a 
model which clearly may be improved upon by perturbation and other methods. 

PACS numbers: 71.23.An,71.20.Dg,71.10.Ca 
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Wigner and Seitz [T][2] originally observed that the major physical properties of sodium 
at one atmosphere (10^ Pa) could be accounted for by considering the crystal as composed 
of two components, namely the ion cores, constituted by the nuclei and their bound elec- 
trons, and itinerant valence electrons, assumed statistically separate from the cores and in 
a paramagnetic state. For the stabilizing volume they found that the cores exerted only 
a weak perturbation on the valence electrons (leading to the nearly free electron (NFE) 
approximation) the core then occupying a relatively small fraction of the unit cell. As a 
starting point for the theory of the metallic state under ordinary conditions, the Wigner 
and Seitz description has been very successful over the years. However, high pressure ex- 
periments conducted recently are now challenging this point of view: when the cores are 
induced to occupy an increasingly larger fraction of the unit cell the indications are that a 
new paradigm may be appropriate, as is suggested here. 

According to system, pressure can reduce linear dimensions by as much as 50% and under 
these conditions the alkalis depart notably from their expected "simple metal" (or NFE) 
behaviors [5| . Their crystal structures at room temperature generally proceed from BCC to 
FCC, and then to non-close packed (for a concise exposition of the structures adopted by the 
alkalis, see [1] and references therein). The latter are difficult to understand intuitively within 
the simple NFE and nuclear- centric viewpoint of Wigner and Seitz. Theoretical work also 
indicates a breakdown of the NFE model: Neaton et al. [5] |6], among others, reported ab 
initio calculations for the band structure of lithium and sodium at high compression. As they 
showed, the bands are far from NFE-like and the occupied band-width is much smaller than 
the familiar SCl/r^ eV expected for free electrons. They also noted |6] that the combined 
effects of Coulomb repulsion, Pauli exclusion, and orthogonality result in an increase of 
valence electron density in the interstitial regions the valence electrons evidently being forced 
away from the near core regions characteristic of nuclear-centric electron distributions. 

This behavior is shown to follow from a simple but quite different physical model which 
focusses on the increasing role played by core exclusion (and near neighbor core proximity) 
upon systematic increase of density, at least at energies typified by Fermi energies. First, 
consider a system of N valence electrons and N spheres (the "ions") occupying a common 
volume Q. Though many of the points to be made are general, the nature of the problem is 
illustrated with two simple lattices, here cubic Bravais lattices (BCC and FCC), with sites 
{R}, and the cores will ultimately be taken to rigorously exclude the valence electrons. As 
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(a)FCC (b)BCC 

FIG. 1: Interstitial Wigner Seitz cells. The red surfaces represent excluding spheres at close 
packing (one sphere is removed in the FCC cell for visibility). Indicated in blue is the plane in 
which the density for Figure |3] is computed and the green cylinder in the FCC cell represents the 
domain for Figure |2} 

in the formulation of the NFE model, electrons are also initially taken as non-interacting, 
though correlation effects will be of later significance. We consider the thermodynamic limit 
oo,Q — > oo,Q/N = 47rrfag/3 and use periodic Von Karman boundary conditions on 
dQ. As a starting point, the ion-electron interaction v is modeled as (in Hartrees): 

r , 



v{r) = VoQir, - 



(1) 



where is the Bohr radius, 6 is the Heaviside function, r^ao is the ionic radius and Vq 
reflects how strongly excluding the potential is. The periodic non-interacting electronic 
Hamiltonian is thus: 
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for which the Schrodinger equation is 



(2) 



(3) 



(a paramagnetic ground state is again initially assumed). As with the NFE approach Bloch's 
theorem fixes the general form of solution of ([s]). Using scaled variables: 

Q R r 



n 



{rsaoY 



X 



V = r^aoV, Vo(rs) 
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the dimensionless form of (Pl) is: 
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For fixed fc = rc/vg, two values of Vq lead to renormalized eigenvalue equations independent 
of r^; Vo = and Vq ^ oo, the latter case related to the Quantum Lorentz Gas with a 
periodic configuration of scatterers ^ . The classical limit will have links to certain billiards 
problems. 

The two-dimensional equivalent yields interesting insight into the problem ahead. Con- 
sider a hexagonal lattice of discs also taken as largely impenetrable to the electrons, as may 
be obtained by considering a (111) lattice plane either for BCC or FCC. At the limiting 
packing fraction in the plane of 7r/2v^, the circular regions touch, and per unit cell there 
remain two involuted triangular regions either of which must eventually contain a single elec- 
tron leading to discrete but highly degenerate and eventually polar ground states. These 
triangular regions are essentially isolated and we know the band-width must vanish at close 
packing in this case. As the radius of the discs approach their close packing value, the tri- 
angular regions are connected by narrow windows and finite band-widths must result. The 
proximity of excluding discs (and their associated boundary conditions) will tend to heavily 
suppress the electronic density at these windows, making it more favorable energetically for 
the density to accumulate at the center of a triangular region. The band-width is expected 
to decay rapidly to zero as the radius of the discs approach their close packing value. Cal- 
culations carried out using the method described next show that this is exactly so. In three 
dimensions, it is not possible to fully isolate an interstitial volume, even at close packing; 
however, we expect the same argument to be valid, ie as the length scale associated with 
the windows becomes much smaller than the length scale associated with the inner region 
of the interstitial Wigner Seitz (IWS) cell, the band-width should decrease rapidly. 

The associated eigenvalue problem can be formulated using a plane wave basis. For set 
{K}, the reciprocal lattice, k a point in the first Brillouin Zone (IBZ), and n a band index, 
the eigenvalue problem proceeds from: 
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FIG. 2: Electronic density along the main body diagonal of the IWS cell of the FCC lattice 
for various values of the ratio Tc/r^. The density is normalized to its maximum value at O for 
comparison purposes. T and O are the tetrahedral and octahedral sites. 

yielding the standard 

(k)K'- (8) 

K 

For finite Vq, we consider an initial finite plane wave expansion. Reciprocal lattice vectors 
K are included in complete stars, all having a radius smaller than a cutoff K\. Sums on the 
IBZ are carried out using the special point technique [8] [9]. 

The formal limit Vq oo must eventually be considered as a boundary condition on the 
electronic orbitals, by demanding that they vanish on the surfaces of the excluding spheres. 
Implementing such a condition completely is quite complex and we use an approximate 
alternative: in the limit of large values of Vq, as is physically clear, the results become 
insensitive to the value of Vq and thus approximate the limit Vq oo. Convergence can 
be assessed by computing the lowest eigenvalue at F (the center of the IBZ), at special 
points on the surface of the IBZ (H,N,P for the BCC lattice, X,W,L,K for the FCC lattice) 
and at points k pointing towards these special points and of magnitude |k| = k^p (the free 
electron Fermi wave vector). The parameters Vq and K\ were chosen to produce reasonable 
convergence with moderate computational effort [lOj. 

From the solution for the states ([6]) we determined the electronic density, centering the 
Wigner Seitz cell around an interstitial site (defining an IWS cell, see Figure [l|). In the 
coordinates of the standard cell the octahedral and the tetrahedral points of the FCC lattice 
have coordinates (|, |, |) and (|, |, |) and the interstitial site of the BCC lattice has coor- 
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(a)FCC: r^jr, = 0.5 
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(b)FCC: r./r, = 0.7 
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(c)BCC: r^jr, = 0.5 (d)BCC: r^/r^ = 0.7 

FIG. 3: Electronic density (from (|3])) in the z=0 plane. The intersection between the IWS cell and 
the z=0 plane is outlined in red. Note that on the approach to close packing, the electronic density 
develops maxima at positions in the IWS cell that are maximally distant from the neighboring 
spheres. 

dinates (0,0, |) (these choices are not, of course, unique). The density has been computed 



on a plane of constant z value centered on the interstitial sites (Figure [3]), and also along 
the body diagonal of the FCC IWS cell, from tetrahedral to octahedral to tetrahedral site 
(Figure [2]). It clearly has a localized character and this provides insight into the different 
behavior seen for the two lattices as r^/r^ approaches close packing. As spheres take up in- 
creasing space in the FCC lattice, the only paths linking cavities centered at the octahedral 
sites of neighboring IWS cells must go through their shared tetrahedral site. As can be seen 
in Figure [2| the boundary conditions efficiently suppress the density at the tetrahedral site 
as TcIts approaches close packing, thus effectively turning the inner surface of the FCC IWS 
cell [it the surface of the Wigner Seitz cell outside the spheres) into a surface where density 
vanishes. Boundary conditions which impose a vanishing of the electronic density at the 



6 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



r /r 

c' s 



FIG. 4: Computed band-width normalized to the free electron Fermi energy as a function of the 
ratio Tc/rs- Inset: absolute band-width, normalized to ec = ^%{rs = Tc). We observe that the 
absolute band-width of the FCC lattice has a turnover around ~ 2rc, whereas for the BCC case 
it does not. 

surface of the cell imply fiat bands. The density at the surface of the IWS cell of the BCC 
lattice, on the other hand, is not as dramatically suppressed even at close packing, allowing 
for a significant band-width. 

Bands and band-widths have been computed in both the FCC and BCC lattices and 
representative results can be seen in Figures |4] and |5j Relative band- widths decrease mono- 
tonically from the free electron value at rd = (which we take to yield the empty-lattice 
bands) to a finite value in the BCC case and to low values in the FCC case at close packing 
{rc/rg — 0.905 for FCC and Tc/ts — 0.879 for BCC) . Furthermore, it appears that the 
Fermi surface develops necks at the L points of the BCC IBZ, and almost so at N in the 
FCC IBZ [11]. In order to further substantiate these results, we also computed the eigenval- 
ues (in units of e^) at F and at H for the BCC lattice using an expansion based on Kubic 
harmonics |12j. We determined the lowest eigenvalue at F in the a representation and the 
lowest eigenvalues at H in the 7 and a representations [13]. The difference between the 7 
eigenvalue at H and the a eigenvalue at F decreases for < Tc/ts < 0.6, suggesting strongly 
that the relative Fermi energy (the band- width) must follow suit. The method produced 
qualitatively different results according to the number of basis functions used at larger ra- 
tios, suggesting that the Bethe Von de Lage method, which tacitly assumes a nuclear-centric 
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(a)FCC: r^/r^ = 0.5 



(b)FCC: r./r, = 0.7 





(c)BCC: rjrs = 0.5 



(d)BCC: r^/rs = 0.7 



FIG. 5: Eigenvalue spectrum for a few values of rc/r^ for excluding spheres in the FCC and BCC 
lattices. Energies are given in units of e^, the free electron Fermi energy, and the energy at the 
bottom of the lowest band is set to zero. The Fermi energy and the free electron Fermi energy (1 
in the graphs' units) are indicated. 

view, does not converge quickly as the ratio approaches close packing. 

The model simulates core exclusion in the simplest way and proceeds from weak coupling, 
NFE like behavior at small relative core volume to narrowing relative band-width (and 
narrowing absolute band-width for FCC) at large relative core volume. The rate at which 
the bands narrow depends on the connectedness of the IWS cell of the lattice and is expected 
therefore to be a general feature. The model reproduces qualitative features observed by 
others [5] [6] , such as the narrowing of the band- width and strong modulations of the electronic 
density leading to interstitial site maxima. It also provides insight into the seemingly strange 
behavior observed earlier and is suggested as a possible initiating paradigm for the alkalis 
under high pressure which may serve as an alternative to the NFE approach. The ionic 
cores are not of course rigorously excluding (as is known from the existence of penetrating 
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states), however the features of the model should partially subsist so long as the average 
kinetic energy of the valence electrons does not exceed the excluding core potential. In that 
sense, it provides insight into ab initio calculations using valence one pseudopotentials with 
repulsive ionic cores. It is unlikely however that the physical system would be well described 
by pseudopotentials designed only to reproduce zero pressure properties when the ionic cores 
are brought to the point of close packing. 

For the physical situation of interest here it becomes clear that the electrons distribute 
themselves in a manner quite different from what the NFE limit would suggest. This can 
be deduced by invoking the simplifying features of the Wigner-Seitz approach, but for the 
largely confining volumes, themselves initially regarded as spheres. As a first approximation 
these have radii = rs{l — (rc/rs)^}^ on the surfaces of which the wave-functions will 
vanish. The corresponding density is proportional to {sin 7r(r/aofs)/7r(r/aofs)}^ and gives, 
near close packing, a quite accurate accounting of the density of Figure [2| which follows from 
wave-functions The corresponding scale of kinetic energy, {ji/fgY /2 Ha, immediately 
suggests the Hubbard physics character of this problem: as the localization energy cost 
becomes large, the energetics of lattice distortions which favor increasing orbital overlap 
and reduced band energy may lead to structural phase transitions. 

Notice also that the corresponding charge distributions will assume even a partially ionic 
form while preserving an essential itinerant or metallic character. Re-introduction of the 
Coulomb interactions is not expected to significantly modify the density profile, however 
the approach will lead to a quite different viewpoint on the character of the effective ion- 
ion interactions. The latter should adopt a somewhat longer ranged form than may be 
expected from linear screening arguments in the NFE limit. This has obvious implications 
for collective excitations and also for ionic dynamics (the physical character of the liquid 
state and melting problem, for example). The behavior of the electronic density also sheds 
light on the structural phase transitions which occur at high pressure: the boundary con- 
ditions in the FCC lattice lead to a depletion of the density at the tetrahedral sites, thus 
reducing the effective available volume for the electrons beyond what is already taken away 
by the impenetrable spheres: a structural phase transition to a state which avoids such 
depleted small enclosures and allows the electron liquid to reduce its zero-point energy by 
fully occuping available space thus may become energetically competitive at high pressure. 
In particular, planar structures allowing the aggregation of the itinerant charge in effectively 



9 



two dimensional NFE gases in planar interstices may become favored. The optical response 
is also anticipated to be unlike that of NFE-based systems, though a Drude edge should 
be preserved. Because of band narrowing, the density of states at the Fermi level could be 
high, and the possibility within a dynamic ionic environment of electron pairing may be 
enhanced. Some of these features as well as mixtures and complex structures are presently 
under examination. 
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